Circumstantial evidence for a non-Maxwellian plasma 
from femtosecond laser-matter interaction 
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We study ion acceleration mechanisms in laser-plasma interactions using neutron spectroscopy. 
We consider diflerent types of ion-collision mechanisms in the plasma, which cause the angular 
anisotropy of the observed neutron spectra. These include the collisions between an ion in the 
plasma and an ion in the target, and the collisions between two ions in the hot plasma. By analyzing 
the proton spectra, we suggest that the laser-generated plasma consists of at least two components, 
one of which collectively accelerated and can also produce anisotropy in the angular distribution of 
fusion neutrons. 
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I. INTRODUCTION 

The development of small-scale high-intensity laser- 
systems with the chirped-pulse amplification (CPA) tech- 
nique has opened new research fields of the laser-plasma 
interaction [1[ . One of the promising applications of these 
studies is the ion-beam generation from laser irradiation 
on solid targets 0, d, It is reported that protons 
in Mylar (H8C10O4) target irradiation are accelerated 
more effectively than in foil target irradiation. Under- 
standing the ion acceleration mechanism in the laser- 
generated plasma is essential for applications. In this 
connection, nuclear reactions induced by laser-irradiation 
give a unique clue in understanding the ion-acceleration 
mechanism [5]. By replacing the protons in the plastic 
CH„ target by deuterons (CD„), a plasma of deuterium 
ions is generated. In the plasma the reaction D((i, n)^He 
with a Q-value of 3.26 MeV is induced [1,0, 9 ^.nd pro- 
duces monochromatic neutrons. The angular distribu- 
tion of the neutrons showspeculiar anisotropy not only 
on CD-plastic target 0, H, P| but also on D2-gas jet [lo| 
and on both D2 and CD2 clusters [111, [i3 • The observed 
neutron angular distribution gives a direct hint to un- 
derstand ion acceleration mechanisms in aneutronic re- 
actions driven by laser as well 13, l3, [Tsj. 

In a recent paper Habara et al. discuss the re- 
sults of an experimental analysis of the neutron spectra 
in nuclear reactions induced by a laser-irradiation on a 
plastic CD target 50 iim thick. They observed that neu- 
tron counts at 1=1 xlO'^^ W/cm^ is larger and more 
anisotropic than that at 1=2 x 10^^ W/cm^. This is at- 
tributed to the fact that higher intensity laser-pulse can 
accelerate ions more efficiently. The ion temperatures 
of 70 and 300 keV at 1=2 xlO^^ and Ix lO^^ W/cm^, 
respectively, with a similar number of accelerated ions 
{Ni = 10^'^) for both intensities are deduced, using a three- 
dimensional Monte Carlo (3D MC) code. They simulate 
ion acceleration processes under different assumptions 
and conclude that the deuterium ions are accelerated into 
the target and cause the nuclear reaction in the target. 
The directionality of the plasma beam is deduced from 



the comparison to the differential cross section data of 
the reaction D(d, n)'^He in conventional laboratory beam- 
target experiments [H, [H [H [HHO, [H . In this paper 
we consider different types of ion-collision mechanisms in 
plasma, which cause the angular anisotropy of the ob- 
served neutron spectra, including collisions between an 
ion in the hot plasma and an ion in the target (HT), 
and collisions between two ions in the hot plasma (HH), 
using at first the total number of accelerated ions and 
the plasma temperature given in Ref. [8|. This assump- 
tion results in the overestimate of the absolute value of 
the neutron yield, if fusion for collisions between ions 
in the hot plasma are properly included. This compo- 
nent was ignored in Ref. Q. Using the SRIM code [22| 
we estimate the HT component which suggests a smaller 
number of plasma ions compared to Ref. _8J. This, in 
turn, reduces the number of fusion originating from the 
collisions among hot ions in the plasma reconciling to the 
experimental observation. However, this is not the only 
possible explanation for the observed angular anisotropy 
in the neutron data. In fact the angular anisotropy in the 
neutron spectra can be observed, if a part of the plasma 
is collectively accelerated and even in the absence of HT 
mechanism. In order to shed some light on this_point 
we study the plasma distribution reported in [sl, |J and 
show that indeed such a collective component is observed. 
We mention that the origin and the mechanism of accel- 
erated ions have been discussed in detail, for a review 
Ref. [i,[2|, and by now it is known the existence of at 
least two types of ion acceleration mechanisms, i.e., from 
the target front side into the target and from the tar- 
get rear side to the vacuum. One of these mechanism 
becomes predominant depending on the target material 
and thickness or laser parameters. The later can be a 
candidate for the collectively accelerated plasma. Finally, 
we stress the importance of knowing both characteristics 
of fusion product, i.e., the spectra of neutron yield and 
the spectra of plasma ions, under common experimen- 
tal conditions. At present those data are available only 
separately. 
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II. NUMBER OF PLASMA IONS DERIVED 
FROM THE OBSERVED NEUTRON YIELD 

In practice we consider the following two types of 
mechanisms for neutron generation in high intensity laser 
irradiation. 

(A) Collisions between two ions in the laser-heated 
plasma. Both ions are moving with thermal ve- 
locity. Under this assumption the direction of 
the incident reaction channel is random, hence the 
angular distribution of reaction products will be 
isotropic 7] . The contribution to the neutron yield 
from this mechanism is called "HH" in this paper. 

(B) Collisions between an accelerated ion in the laser- 
produced plasma and a cold nucleus in the bulk 
of the target. Under this second assumption the 
angular distribution of reaction products is possibly 
anisotropic. The contribution from this mechanism 
is called "HT" component. 

We stress the importance of considering both mechanisms 
mentioned above, comprehensively, because either mech- 
anism might be predominant, depending on the char- 
acteristics of the target. As an example in the case of 
neutron yield observation from laser pulses irradiation 
on deuterated clusters [11] both mechanisms play a key 
role. 

If we assume that neutrons are produced by the colli- 
sion of the ions in the hot plasma component (HH), in 
terms of the number of the accelerated ions Nj, the num- 
ber of fusion per solid angle, or reaction rate [2J|, is given 
by 

M""^ 1 f 

= —N,n„T / a{v)v^{v)dv^ , (1) 

where tier =10^^/A^ is the plasma critical density [25l.[26t: 
r is the laser pulse duration; u(v) and v are the reaction 
cross section and the relative velocity of the colliding ions. 
In general the reaction cross section <t{v) is given as a 
function of the incident energy, instead of the velocity, 
but here we have written it as a function of the velocity 
to keep the consistency in the velocity integral. Later a 
will be represented as a function of the incident energy. 
We stress that assuming the critical density, which is the 
lowest limit to the real density reached in the experiment, 
the fusion yield given by Eq.([T]) is underestimated. The 
density profile of the plasma simulated in Ref 8], using 
the PIC code, has an exponential shape which varies from 
Aricr to O.lricr- 4'{'v) is the relative velocity spectrum of 
a pair of ions and is given by a Maxwellian-distribution 
at the temperature kTnH =70 or 300 keV: 

where /i is the reduced mass of ions. Eq. ([1]) gives 1.x 
10^ and 4.x 10^, per solid angle, at the temperature of 70 



and of 300 keV, respectively, see Tab]!] The estimated 
yield is comparable with the neutron spectra in Fig.s 3 
and 5 in Ref. Q. The contribution from HH component 
is, therefore, expected to be seen in the figure as a peak 
at the neutron energy 2.45 MeV and to be isotropic 
In order to compare this to their energy distribution, we 
can roughly assume that the neutron distribution is a 
Gaussian distribution with a center-of-mass (CM) energy 
at 2.45 MeV and a width given by the temperature of 
the plasma. This gives an estimate at all angles of 1.4 
xlO^ and 1.3 x 10^ ion/MeV/sr respectively, which is 
seen neither in figures 3 nor 5 in Ref. [11 . In the figures, 
if anything, one sees a clear angular dependence of the 
neutron yield and shifts of the observed peaks from the 
expected energy 2.45 MeV. This implies that either their 
estimated temperature or the number of accelerated ions 
given is too large. In other words the authors of Ref. Q 
should have estimated the number of neutrons coming 
from the HH component, and show that this component 
is negligible compared to the HT contribution. As we 
will show in the following discussion, the contribution 
from the HT component, which has the correct angular 
dependence of the observed neutron spectra, is indeed 
dominant but with a smaller number of plasma ions. 

If we consider the collisions between the ions in the 
plasma and the almost stable nuclei in the target (HT), 
the angular distribution of reaction products is expected 
to be anisotropic. To estimate the neutron yield from the 
HT component, one should take into account that one 
of the colliding ions is at rest in the laboratory frame. 
Therefore in the reaction rate per pair of colliding ions, 
the velocity spectrum Eq.([2]) is modified as: 

where mi is the mass of ions in plasma, instead of the 
reduced mass. One can define an effective temperature 

as, 

fcr^T = {t^/mi)kTHH- (4) 

Now we can use the effective temperature to estimate 
the number of fusion. For simplicity we estimate the 
most probable energy of the plasma ions that cause the 
nuclear reaction given by the Gamow peak energy {Eq) 
[2^ . The Gamow energy can be found using the saddle 
point method, i.e.: 

where b — 3I.28Z1Z2M5 (keV^), denoting the atomic 
number of the colliding nuclei Zi , Z2 and the reduced 
mass number M. We remind that M — AiA2/{Ai+ A2), 
where Ai and A2 are the mass numbers of the colliding 
ions, respectively, is different from the reduced mass /i. 
The temperature kT is replaced by the plasma tempera- 
ture, i.e., the temperature of the HH component in con- 
ventional discussions, but for the present case of the HT 
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component, kT is replaced by the effective temperature. 
Then, 
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Eg = 



(6) 



We speciahy mention that Eq. ([5|) is vahd in the case of 
sub-barrier reactions, because this condition comes from 
the product of a MaxweUian and the Coulomb barrier 
penetrability [1^ . The height of the Coulomb barrier for 
the reaction D(d, n)^He is estimated to be about 470 keV. 
By approximating the energy of the accelerated ions by 
the Gamow energy Eq. ([6]), the neutron yield per solid 
angle is written as: 
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where cr{EG) is the reaction cross section at the Gamow 
energy. The reaction cross section data at this energy 
are taken from NACRE compilation [27*] and given in 
Table HI The number density of the solid polyethylene 
target: tlt = 1.2 x 10^'^ atoms/cm'^ and d is the pro- 
jected range [l^l of the accelerated ions in the target. 
Using the above listed numbers the neutron yield per 
solid angle is also reported in Tab. [J As a consequence 
of this simple estimate, the yield per solid angle under 
this assumption is a factor 10^ higher than the yield from 
the collisions between two ions in the plasma. This is 
because the number density of the solid target is much 
higher than that of the laser-generated plasma. However 
from eq. (jl]) the effective HT temperature is lower than 
the HH one, depending on the asymmetry in mass of the 
fusion ions. The higher temperature in the HH collision 
mechanism has an advantage of the higher fusion cross 
section. There is, therefore, a competition between two 
mechanism depending on the laser-intensity. This fea- 
ture might explain why at lower intensities (i.e. lower T) 
the neutron angular distribution is less anisotropic: At 
the lower intensity irradiation, the contribution from the 
HT component to the neutron yield is suppressed com- 
pared with the higher intensity irradiation. In fact the 
temperature of the HT component is low and most re- 
actions happen below the Coulomb barrier. The nuclear 
reactions below the barrier are exponentially suppressed. 
Therefore the HH mechanism, which has higher temper- 
ature, contributes to the less anisotropic neutron yield 
angular distribution. On the other hand, if the plasma 
temperature is high enough, the corresponding Gamow 
energy is above the Coulomb barrier. At energies above 
the Coulomb barrier, i.e. high T, to have higher densities 
gives higher fusion probabilities (above the barrier fusion 
probabilities depend quadratically on densities and do 
not depend much on T) thus collisions between an ion 
from the plasma and an ion from the target becomes 
more probable. We expect that in some temperature re- 
gion there is a transition from collisions occurring mainly 
in the plasma to collisions occurring between the plasma 
ions and the target ones. An experimental detailed inves- 
tigation of this transition region would be very interesting 
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FIG. 1: Experimental data of proton energy spectra for Mylar 
targets 6 (asterisks) and 13 (crosses) /xm-thick retrieved from 
Lee et al. [4|. The thick curves represent fitting of the proton 
spectra using MaxweUian distributions and distributions with 
extra accelerations. The thin curves are the different compo- 
nents and are summed into the thick curves. The obtained 
fitting parameters are given in Tab. |IT1 



and instructive in order to understand the microscopic 
dynamics of fusion. 

The HT contribution to the neutron yield is a factor 
10^ higher than experimentally observed data Q, as well. 
This again implies that either the estimated tempera- 
tures are too high or the estimated number of ions in the 
plasma is too large. From this simple estimate and the 
one above from the HH component we can argue that the 
number of ions in the plasma could be at least two orders 
of magnitude smaller than that estimated in Ref. Q, i.e., 
about 10^^, instead of 10^^. This reduction of the number 
of accelerated ions results in a suppression of the neutron 
yield from the HH component which becomes negligible. 
A more involved calculation solving eq.(l) for the HT 
component numerically gives results in agreement with 
our simple estimate and will be discussed in more detail 
in a following paper [2^. 



III. PLASMA TEMPERATURE SUGGESTED 
BY THE PROTON SPECTRA FROM PLASTIC 
TARGETS IRRADIATION 

Another evidence of the smaller number of plasma ions 
is shown by a direct observation of proton spectra in laser 
irradiation Fig. 1 in Lee et al. [i^ shows the proton 

energy distribution at the laser intensity / — 2.2x 10^^ 
W/cm^ similar to Q, and the target material is Mylar 
or aluminum. We are especially interested in the results 
from the Mylar target, because the characteristics of the 
produced plasma should be close to the characteristics 
of the plasma from the deuterated plastic targets. We, 
therefore, selected two results from 6 and 13 ^m thick 
Mylar target. A peculiar feature of the proton energy 
distribution is that it exhibits bumps in the higher energy 
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TABLE I: Plasma temperature (kTHn) and neutron yield per solid angle {Nj: / Sr) from the HH component and the effective 
temperature (fcT^'t/) and Nj,^'^'' /Sr from the HT component at the given laser intensities. 



I (W/cm^) 


kTHH (keV) 




fcr^^T^'(keV) 


£'G(keV) 


(j(£'G)(10"^''cm2) 


d(/im) 




2x10'*" 


70 




35 


67 


24 


0.8 


4. X 10'' 


1x10^*^ 


300 


4.x 10'' 


150 


176 


66 


2. 


3.x 10" 



TABLE IL Fitting parameters for proton energy spectra from 
6 and 13 nm thick targets irradiation. Especially for 13 fj,m 
thick target we have used three components, two of which 
have different relative velocities with respect to the hot- 
plasma component. 



T thick. 




kTHH 


Cl 


kTi 


Eo 


(/im) 


(10") 


(keV) 




(keV) 


(keV) 


6 


(3.0±0.2) 


44.±1. 


(3.5±0.7) lO'' 


52.±5. 


170. ±20. 


13 


(3.5±0.1) 


61. ±1. 


(3. ± 2.) 10^ 


78.±30. 


210. ±60 








(4. ± 1.) 10^ 


14.±5. 


695.±10 



region. This feature could be attributed to the existence 
of at least two different components in the plasma. To 
describe these characteristics, we use a Maxwellian dis- 
tribution in terms of the energy of the protons in the lab. 
system {E): 
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To this we add a Maxwellian distribution of a moving 
source [291: 



d^N Cl VeW 
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where E' — E — 2\/EEa + Eq. Both equations are nor- 
malized to l/47r. The latter gives a part of the plasma in 
translational motion with respect to the laboratory sys- 
tem with a collective kinetic energy Eq. The obtained 
fitting curves are shown by thick solid and dashed curves 
for the spectra at 6 and 13 /im thick targets irradiation, 
respectively, in Fig. [TJ The thin curves show each com- 
ponent which add into the thick curves. Especially for 
the spectra at 13 ^m thick target, we have used a sum 
of three components, two of them have different extra 
acceleration energies Eq. The corresponding fitting pa- 
rameters are summarized in Tab. HT] From the fitting 
we can deduce the number of ions in the plasma iV^ ~ 
3x10^^ and the plasma temperature of UThh ~ 44 and 
61 keV for 6 and 13 /im thick targets, respectively. The 
deduced number of plasma ions is in agreement with our 
yield estimation from the neutron measurement. More- 
over one can see clearly that the proton spectra differ 
from a Maxwellian but a part of plasma is accelerated 
to higher energy with Eq. Similar features of the pro- 
ton spectra in laser-produced plasma can be found in 
Spencer et al. as well. The possible reasons for such 



a collective acceleration will be discussed in a following 
publication j28i] . 

In Fig. 7 and 8 of Q, the authors discuss the ratio of 
the neutron yields at the angle of 23° to that at the angle 
of 67° from the target rear normal. They found that the 
ratio can be larger than 3. In Fig. [5] the ratio of the 
differential cross sections in beam-target experiments [l^, 
at laboratory angles 23° to 67° is shown 
by thick curves. There is a discrepancy between [sc72] 
data [iSj] and [1173] [l^l in the energy region higher than 
3 MeV. Nevertheless the figure shows clearly that the 
ratio can, indeed, be larger than 3 and less than 3.5 at the 
incident deuteron energy El from 1.75 to 1.96 MeV. Thus 
the accelerated deuterons which contribute to the nuclear 
reactions have energies slightly less than 2 MeV. In the 
CM frame E ^0.9 MeV, i.e., the corresponding Gamow 
peak energy is about 900 keV. If we use i?G=900 keV, 
the neutron yield per solid angle from the HT collision is 
estimated to be of the order of 10'' from Eq. ([7]), adopting 
the number of accelerated ions given in This number 
is 10 times higher than that in Table [H In other words, 
to reproduce the experimental data, one needs to assume 
Ni ^^10^°, as the number of accelerated ions in plasma. 
In passing we note that the ratio of the neutron yield 
taken at the angle of 20° to that of 85° is higher than 
the previous one as shown in Fig. [2] by thin curves with 
smaller points. A careful experimental determination of 
this feature would be very useful. 

We can also derive the plasma temperature corre- 
sponding to £^G=900 keV. Since at this energy we are 
above the Coulomb barrier, we can estimate it from the 
classical relation: 



Eg - {■i/2)kTHH, 



(10) 



which gives fcrff/f=600 keV, i.e., higher than the es- 
timated temperature in Ref. This contradiction 
might be solved by considering that a part of the plasma 
deuterons is collectively accelerated at energy Bq, as we 
have shown in the analysis of proton spectra. This ex- 
tra acceleration energy can reach about 700 keV which is 
close to 900 keV. In the presence of an extra acceleration, 
the plasma temperature cannot be derived simply from 
the relation (|10p . The difference between the plasma tem- 
perature deduced from the ratio at two angles and the 
one estimated by Ref. Q suggests that the plasma spec- 
tra is different from a usual Maxwellian distribution. In 
other words, the difference justifies the presence of the 
collective motion of a part of the plasma. Indeed the col- 



5 



c 
o 



3 



1^ 

CD 



CO 
CM 



O 
03 



^ 



th66 — o--- 

sc72 — ■— 

ja77 -H3— 

hu49{neutron) • 

Ii73 ■■■■><■■- 




0.01 



0.1 1 
El (MeV) 



10 



FIG. 2: The ratio of differential cross sections at two angles. 
The ratio at angles 23° to 67° is shown by large points con- 
nected by curves, while the ratio at angles 20° to 85° is shown 
by smaller points connected by thin curves, for a comparison. 
The experimental data are retrieved from EXFOR-data sys- 
tem |30|. hu49-data (pink filled circles) are evaluated from 
the angular distribution of the DD neutron yield [igI ] and 
the otherSjth66-data (grey open circles) [ITII . sc72 (red filled 
squares) [3 1 ja77 (green open squares) [2^ and 1173 (blue 
crosses) [ly], are converted from the "^He angular distribution 
of the DD reaction. 



lisions between two ions in the plasma which is moving at 
a collective energy Eq in the laboratory frame can result 



in the angular anisotropy of the neutron yield. This is 
a possible mechanism to explain the angular anisotropy 
experimentally observed. 



IV. CONCLUSIONS 

In conclusion, we have discussed different possible ion- 
collision mechanisms, which can result in the observed 
anisotropic neutron spectra. When analyzing the proton 
spectra, we have found that there are at least two plasma 
components: one is approximated by a Maxwellian dis- 
tribution and the other has a collective motion relative 
to the former component. Comparing the ratios in the 
neutron counts at two angles in the experiment to the dif- 
ferential cross sections measured in conventional beam- 
target experiments, the most effective energy is estimated 
to be 0.9 MeV with corresponding plasma temperature 
of 600 keV, which is higher than the estimated plasma 
temperature in Ref . [§| . We have discussed a possible so- 
lution of this contradiction, in connection with the collec- 
tive motion of a part of the plasma, which can in princi- 
ple explain the observed neutron angular anisotropy. We 
suggest that an experimental determination of the neu- 
tron angular anisotropy in coincidence with the plasma 
distribution would be very interesting and give useful in- 
formation on the mechanisms at play. 
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